Dynamic properties of atomic collective decay in cavity quantum electrodynamics
Han Yu-Feng1, 2, Zhu Cheng-Jie1, †, Huang Xian-Shan2, Yang Ya-Ping1, ‡
Key Laboratory of Advanced Micro-Structured Materials of the Ministry of Education, School of Physics Science and Engineering, Tongji University, Shanghai 200092, China
School of Mathematics and Physics, Anhui University of Technology, Ma’anshan 243032, China

 

† Corresponding author. E-mail: 1510574@tongji.edu.cn yang_yaping@tongji.edu.cn

Project supported by the National Natural Science Foundation of China (Grant Nos. 11504272, 11774262, 11474003, and 11504003), the National Key Basic Research Special Foundation (Grant No. 2016YFA0302800), the Joint Fund of the National Natural Science Foundation of China (Grant No. U1330203), the Fund from the Shanghai Science and Technology Committee (STCSM) (Grant No. 18JC1410900), and the Natural Science Foundation of Anhui Province, China (Grant Nos. 1408085MA19 and 1608085ME102).

Abstract

We theoretically study the collective decay of two atoms trapped in a single mode cavity and we describe the evolution of the population of Dicke states. We show that the collective decay property is strongly dependent on the phase of atomic radiation and the speeding up of collective decay can only be observed in a bad cavity regime. For in- or out-phase case, this occurs due to the quantum interference enhancement, no matter which atom is excited initially. For π/2 phase, the speeding up of collective decay takes place if the first atom is excited at the beginning. However, it disappears due to the quantum interference cancellation if the second atom is excited. Compared with the in-phase and out-phase cases, we also show that the speeding up of collective decay can be significantly enhanced in strong coupling regime for π/2 phase, although one atom is decoupled to the cavity in this condition. The study presented here is helpful to understand the physical mechanism of collective decay in cavity quantum electrodynamics and it provides a useful method to control the collective decay phenomenon via quantum interference effect.

1. Introduction

Quantum interference effects are basic features of quantum mechanics, they are also the foundation of quantum optics. There are various interesting optical phenomena that are essentially based on quantum interference effects, for example: electromagnetically induced transparency,[1] lasing without inversion (LWI),[2,3] and refractive index enhancement without absorption.[4] In addition, quantum interference can lead to the modification of atomic spontaneous emission, such as, spontaneous emission enhancement in photonic crystals[5,6] and reduction even cancellation in three-level atoms.[79]

It is well-known that the spontaneous emission rate emitted by quantum emitters can be modified by inhomogeneity of the cavity environment in the Purcell effect.[10] Purcell also demonstrated that the spontaneous emission rate is not an intrinsic property of the emitter. Since then, great attention has been drawn on resonant systems to study the enhancement of the spontaneous emission decay rate, such as plasmonic waveguides,[11] photonic crystals,[6,12,13] nanocavities,[14,15] and nanoshells.[16] However, it is hard to enhance the total (collective) emission in an ensemble of quantum emitters system due to the severe condition of the emitters’ accurate location as well as the neighboring large and homogeneous electric field distributions in resonating systems.

Collective emission or decay occurs in an ensemble of quantum emitters system where the emitters are coupled to a lossy cavity[17,18] or close to each other.[19,20] Collective decay is of fundamental interest in quantum optics and has attracted tremendous attention[21] since the seminal work by Dicke,[22] especially relating to the phenomenon to super-radiance.[19,2326] The phenomenon can be attributed to the constructive quantum interference between emitted waves and has been widely investigated in various photonic environments.[2733] Through coherent control of collective decay, the super-radiant states can be converted into subradiant ones and vice versa in an extended atomic ensemble.[34] Chen et al. proposed a scheme to realize entanglement swapping via collective decay, which is a procedure to create entanglement between two qubits.[35] Agarwal et al. showed determinstic generation of Werner states of the collective decay dynamics of a pair of neutral atoms coupled to a leaky cavity and strong coherent drive.[36] Wang et al. also presented a simple scheme to effectively generate maximal pure-state entanglement between non-interacting atoms through purely collective decay.[37] Despite all of these research achievements, the nature or physical mechanism of collective decay is, even today, not well understood.

In this paper, we have carefully studied the dynamics of atomic collective decay in cavity quantum electrodynamics. We assume that two atoms are trapped in a single mode cavity with different positions and that one of them is initially excited by a single photon. It is found that the decay properties of the Dicke states are strongly dependent on the phase between two atoms. We show that the speeding up of collective decay can only be observed in the bad cavity regime due to the quantum interference effect. In the case of in-phase or out-phase radiation, the decay rate of the Dicke state is much faster than that of the single atom, which is recognized as a witness of the single photon super-radiance.[38] Here, it should be emphasized that the obtained single photon super-radiance is referred to the collective decay enhancement of two atoms trapped in a cavity via single photon excitation comparing with the spontaneous decay of one atom in vacuum. We also show that the properties of the collective radiation change significantly if the atoms feel different coupling strengths to the cavity. For example, in the case of π/2 phase, the speeding up of the collective decay can be observed when the first atom is excited initially, while it disappears if the second atom is excited at the beginning. Particularly, in the case of π/2 phase, the speeding up of the collective decay can be significantly enhanced by increasing the coupling strength. The results presented in this work are helpful for understanding the nature of the collective decay phenomenon in the cavity and can be implemented in other systems.

2. Bad cavity regime

To begin with, we consider that two two-level atoms with the same resonant transition frequency ωA are weakly coupled to a single mode cavity with frequency ωC [see Fig. 1(a)]. State |g〉 (|e〉) corresponds to the ground (excited) state of each atom, and γ denotes the spontaneous decay rate from state |e〉 to state |g〉. Since the two atoms are located at different positions, the coupling strengths between atoms and cavity are different and denoted by g1 and g2 respectively.

Fig. 1. (color online) Physical model sketch and Dicke state diagram. (a) Two atoms are trapped in a single mode cavity with different positions. The ground (excited) state of each atom is labeled as |g〉 (|e〉). γ and κ represent the decays of atoms and cavity, respectively. The coupling strengths between atom and cavity are denoted by g1 and g2 respectively. (b)–(d) The coupling and decay between the lowest four Dicke states for phases between two atoms ϕz = 0, π/2, and π, respectively. Here, we define the collective states .

The dynamical evolution of this atom–cavity coupled system can be well described by the Jaynes–Cummings Hamiltonian in the frame rotating with the cavity frequency,[39,40] i.e., where ΔA = ωAωC is the detuning of the atomic resonance frequency from the cavity frequency. is the raising (lowering) operator for the j-th two-level atom, and a (a) is the cavity photon creation (annihilation) operator. For simplicity, we can assume that the first atom is localized at the antinode of the standing wave in the cavity, so that the coupling strength between the first atom and the cavity is g1 = g0. Then, the coupling strength between the second atom and the cavity is g2 = g0cos(ϕz) with ϕz = 2πΔz/λC being the phase shift between the radiation of the two atoms. Here, Δz is the distance between two atoms and λC = c0/ωC is the wavelength of the cavity mode where c0 is the light speed in vacuum.

In the normalized unit with ħ = 1, the evolution of the density matrix ρ of this system, including all the dephasing terms, is given by where and originate from the spontaneous decay of each two-level atom at rate 2γ and the decay of the cavity field at rate 2κ. It is obvious that the collective decay rate is phase-dependent and results in the super-radiance and subradiance effects.

In general, equation (2) can only be solved numerically since high-order photon states must be considered in the strong coupling regime. However, for a bad cavity and weak coupling strengths, only few photons in the cavity can interact with the atoms so that the master equation can be simplified by neglecting the cavity component, yielding where are the effective operators for Dicke states.

Considering the lowest three states of this simplified system for single photon excitation, i.e., |1〉 = |gg〉, |2〉 = |eg〉, and |3〉 = |ge〉, we can obtain a set of equations from Eq. (3), which reads where , and for conversation. Here, we must point out that the population of state |ee〉 is zero since there is only one photon in the cavity.

If one assumes with τ being the time scale of interaction, equation (4c.) can be solved under the steady state approximation, i.e.,

Then, by inserting the above equation into Eqs. (4a)–(4b) and using the initial condition ρ33(0) = 1 (i.e., the system is initially prepared in state |3〉 = |ge〉), we can obtain the populations in states |eg〉 and |ge〉, i.e., where with , and .

Equations (6a) and (6b) predict many interesting characteristics for the dynamical evolution of populations in states |eg〉 and |ge〉. For example, it is noted that if the phase between the two atoms satisfy ϕz = π/2, the first atom will not be excited [i.e., Peg(t) = 0], and the population in state |ge〉 decays at rate 2γ [i.e., Pge(t) = exp(−2γt)]. However, in the case of ϕz = 0 or π, we can obtain . Therefore, it is easy to find that Peg(t) ≈ Γt and Pge(t) = exp(Λ+t), which implies that the population in state |eg〉 will experience a linear gain at the beginning, and an enhanced damping of the population in state |ge〉 can be observed.

The physical mechanism can be understood by rewriting the Hamiltonian by using the collective states , and |ee〉 as basis.[41] To this end, the Hamiltonian can be expressed in terms of the collective operators , yielding HI = H+ + H with and g± = g0(1 ± cosϕz). Using the initial condition, we can obtain that at the beginning. For the case of ϕz = 0 (ϕz = π), only |+,0〉 ↔ |gg,1〉 (|-,0〉 ↔ |gg,1〉) coupling is allowed but the |-,0〉 ↔ |gg,1〉 (|+,0〉 ↔ |gg,1〉) coupling is forbidden [see Figs. 1(b) and (d)]. Thus, the single photon super-radiance can be observed via the constructive interference between |+,0〉 ↔ |gg,1〉 → |gg,0〉 (|-,0〉 ↔ |gg,1〉 → |gg,0〉), and |+,0〉 → |gg,0〉 (|-,0〉 → |gg,0〉) pathways. However, for the case of ϕz = π/2, both |+,0〉 ↔ |gg,1〉 and |-,0〉 ↔ |gg,1〉 couplings are allowed [see Fig. 1(c)], but the coupling strengths have opposite sign, which results in an destructive interference and the decay rate of the state |ge〉 remains that of single atom (i.e., 2γ).

To verify this analysis, we numerically solve Eqs. (2) and plot the populations in states |eg〉 [panel (a)] and |ge〉 [panel (b)] as functions of the phase ϕz and the evolution time t in Fig. 2.

Fig. 2. (color online) The populations in states |eg〉 [panel (a)] and |ge〉 [panel (b)] with initial condition Peg(0) = 0 and Pge(0) = 1. The system parameters are chosen as g0 = 1, κ = 3, and ΔA = 0.

The initial conditions are chosen as Peg(0) = 0 and Pge(0) = 1. Here, we choose the system parameters as g0 = 1, κ = 3, and ΔA = 0. If one chooses ϕz = π/2 where the population in state |ge〉 exponentially decays at the rate of 2γ [see panel (b) of Fig. 2], which is the same as the spontaneous emission rate for the single atom. The first atom is always in state |g〉, i.e., Peg(0) = 0 [see panel (a) of Fig. 2]. However, if ϕzπ/2, the population in state |ge〉 drops more quickly than the rate of 2γ and reaches the maximum when ϕz = 0 and π, as shown in panel (b) of 2. Simultaneously, the first atom absorbs the photon emitted from the second atom and jumps to the excited state, i.e., Peg(0) ≠ 0. In the case of in-phase or out-phase radiation, the possibility in the excited state for the first atom can be maximumly achieved due to the constructive interference effect [see panel (a) of Fig. 2].

Next, we will study the collective decay effects if the two atoms are initially prepared in state |eg〉 (i.e., Peg(0) = 1 and Pge(0) = 0). Using the same method, the solutions of Eqs. (4a)–(4c) can be obtained easily, which reads

The solutions are similar to those given in Eqs. (6a) and (6b), but hold different dynamical properties. Assuming ϕz = 0 or π, we can obtain , which implies that the population in state |eg〉 decays at the rate of . In fact, the factor represents the Purcell effect. This enhancement of the decay can also be explained by constructive quantum interference resulting from the additional pathway of |+,0〉 ↔ |gg,1〉 → |gg,0〉 (|-,0〉 ↔ |gg,1〉 → |gg,0〉) [see Figs. 1(b) and 1(d)]. In the case of ϕz = π/2, the population in state |eg〉 is also . Contrary to the case of Pge = 1 as the initial condition, the decay rates of the state |eg〉 are always enhanced, even if only single atom is coupled to the cavity. This feature can be attributed to the constructive quantum interference between the pathways of |+,0〉 ↔ |gg,1〉 and |-,0〉 ↔ |gg,1〉 [see Fig. 1(c)] which have the same coupling strength when ϕz = π/2.

By once more numerically solving Eqs. (2) with initial conditions of Peg(0) = 1 and Pge(0) = 0, we show the dynamical evolution of the populations in states |eg〉 and |ge〉 versus the normalized phase ϕz and the evolution time t in Fig. 3. The system parameters are the same as those used in Fig. 2. It is clear that the decay rate of the state |eg〉 is always enhanced (i.e., bigger than 2γ, which is the spontaneous decay rate in vacuum) and almost independent of the phase between two atoms [see panel (a) of Fig. 3]. In the case of ϕz = π/2, the second atom is always in ground state for the reason of uncoupling with the cavity. However, if ϕzπ/2, the second atom exchanges energy with the first atom via coupling with the cavity and the possibility of detecting the states |ge〉 reaches the maximum in case of ϕz = 0 or ϕz = π [see panel (b) of Fig. 3].

Fig. 3. (color online) The populations in states |eg〉 [panel (a)] and |ge〉 [panel (b)] with initial condition Peg(0) = 1 and Pge(0) = 0. The system parameters are the same as those used in Fig. 2.
3. Bad cavity and strong coupling regime

Now, we consider the case of g0κ,γ, i.e., the strong coupling regime. Since the back reaction becomes important in this case, equation (2) must be solved numerically to study the dynamical evolution of the decay in Dicke states, nevertheless, equation (3) is not effective. To make the results more explicit, only the cases of ϕz = 0 and ϕz = π/2 will be adopted under consideration in the following discussion (the case of ϕz = π is the same as ϕz = 0).

As shown in Fig. 4, we show the dynamical evolution of populations in states |eg〉 and |ge〉 with the coupling strength g0 = 0.01 (the blue curve), g0 = 1 (the green curve), and g0 = 10 (the red curve), respectively. Here, we first take the initial condition as Peg(0) = 0 and Pge(0) = 1. In the case of ϕz = 0, we find that when the coupling strength between the atom and cavity is weak (e.g. g0 = 0.01), the population in the state |ge〉 decays slowly at the approximate rate of 2γ [see blue curves in panel (b)], which is the same as the spontaneous emission rate for the single atom. The first atom cannot absorb the photon emitted from the second atom and remain in the ground state due to weak coupling strength [see the blue curve in panel (a)]. As the coupling strength increases, we show that the population in the state |ge〉 decays far more quickly than spontaneous emission for the single atom in vacuum [see the green curve in panel (b)]. Meanwhile, the possibility of being excited for the first atom by absorbing the photon emitted from the second atom gradually increases [see the green curve in panel (a)]. For strong coupling strength (here we choose g0 = 10), the decay rate of the population in the state |ge〉 increases further. Resulting from the strong coupling, resonance between |+,0〉 ↔ |gg,1〉 makes both Peg and Pge oscillate intensively at the frequency [see the red curves in panels (a) and (b) in Fig. 4]. Moreover, two atoms continuously exchange energy from each other via the cavity. Therefore, the possibility of detecting the states |eg〉 exceeds 50 percent.

Fig. 4. (color online) The populations in states |eg〉 [panels (a) and (c)] and |ge〉 [panels (b) and (d)] with initial conditions of Peg(0) = 0 and Pge(0) = 1. Here, ϕz = 0 for the top two panels, but ϕz = π/2 for the button ones. The coupling strength is chosen as g0 = 0.01 (blue curve), g0 = 1 (green curve), and g0 = 10 (red solid curve). Other system parameters are the same as those in Fig. 2.

In the case of ϕz = π/2, we show that the decay rate of state Pge is the same as that for the single atom, and Peg = 0 no matter how strong the coupling strength is [see panels (c) and (d) in Fig. 4]. This is because the quantum destructive interference takes place between the |+,0〉 ↔ |gg,1〉 and |-,0〉 ↔ |gg,1〉 couplings. We must point out that only in this case is the second atom transparent to the cavity.

In Fig. 5, we choose Peg(0) = 1 and Pge(0) = 0 as the initial condition. Other system parameters are the same as those used in Fig. 4. It is clear that, in the case of ϕz = 0, the decay rate for Peg increases with the increasing coupling strength. For strong coupling, Peg demonstrates damped oscillation with the frequency of [see panel (a)]. At the same time, we discover that the possibility of being excited for the second atom by absorbing the photon emitted from the first atom also increases with the increasing coupling strength [see the green curve in panel (b)]. For strong coupling, Pge also exhibits damped oscillation with the same frequency with Peg [see the red curve in panel (b)]. Now, the decay rate for the first atom increases with the increasing coupling strength, which is the well-known Purcell effect.[9] In the case of ϕz = π/2, Peg shows Rabi oscillation but with the frequency of 2π/(2g0) [see panel (c)]. However, the population of the state |ge〉 is always zero [see panel (d)] due to destructive interference effect.

Fig. 5. (color online) The populations in states |eg〉 [panels (a) and (c)] and |ge〉 [panels (b) and (d)] with initial conditions of Peg(0) = 1 and Pge(0) = 0. Here, ϕz = 0 for the top two panels, but ϕz = π/2 for the button ones. The coupling strength is chosen as g0 = 0.01 (blue curve), g0 = 1 (green curve), and g0 = 10 (red solid curve). Other system parameters are the same as those in Fig. 2.
4. Conclusion

We have shown that the behavior of the atomic collective decay significantly changes with the phase between the two atoms trapped in a high damping cavity. This arises from the quantum interference between Dicke states. We found that in the case of in-phase or out-phase radiation, the speeding up of the collective decay occurs due to the constructive quantum interference, no matter which atom is excited initially. In the case of π/2 phase, the speeding up of the collective decay can also be observed if the first atom is excited initially. Although in this case the second atom does not couple to the cavity, the constructive interference still exists and results in an enhanced decay rate of the Dicke state. However, the speeding up of the collective decay disappears due to the destructive quantum interference if the second atom is excited at the beginning. In addition, the speeding up of the collective decay can be significantly enhanced in the strong coupling regime for the π/2 phase. The study presented here is helpful to understand the nature of the speeding up of the collective decay in cavity quantum electrodynamics, and provide an effective method to control the collective decay via quantum interference effect.

Reference
[1] Harris S E 1997 Phys. Today 50 36
[2] Harris S E 1989 Phys. Rev. Lett. 62 1033
[3] Scully M O Zhu S Y Gavrielides A 1989 Phys. Rev. Lett. 62 2813
[4] Fleischhauer M Keitel C H Scully M O Su C Ulrich B Zhu S Y 1992 Phys. Rev. 46 1468
[5] Yang Y Zhu S Y 2000 Phys. Rev. 61 043809
[6] Liu Z Q Feng T H Dai Q F Wu L J Lan S Ding C R Wang H Z 2010 Chin. Phys. 19 114210
[7] Scully M O Zubairy M S 1988 Opt. Commun. 66 303
[8] Ansari N A Gea-Banacloche J Zubairy M S 1990 Phys. Rev. 41 5179
[9] Zhu S Y Chan R C Lee C P 1995 Phys. Rev. 52 710
[10] Purcell E M Torrey H C Pound R V 1946 Phys. Rev. 69 37
[11] Sokhoyan R Atwater H A 2013 Opt. Express 21 32279
[12] Noda S Fujita Asano M T 2007 Nat. Photon. 1 449
[13] Li J Wan R G Yao Z H 2016 Chin. Phys. 25 104204
[14] Russell K J Liu T L Cui S Hu E L 2012 Nat. Photon. 6 459
[15] Qiu L Zhang K Li Z Y 2013 Chin. Phys. 22 094207
[16] Tam F Goodrich G P Johnson B R Halas N J 2007 Nano Lett. 7 496
[17] Bohnet J G Chen Z Weiner J M Meiser D Holland M J Thompson J K 2012 Nature 484 78
[18] Bonifacio R Schwendimann P Haake F 1971 Phys. Rev. 4 302
[19] Gross M Haroche S 1982 Phys. Rep. 93 301
[20] Lin G D Yelin S F 2012 Advances in Atomic, Molecular, and Optical Physics Netherlands Elsevier 295 329 10.1016/B978-0-12-396482-3.00006-5
[21] Hesham F Enrica C 2017 Chin. Phys. Lett. 34 114101
[22] Dicke R H 1954 Phys. Rev. 93 99
[23] Men'shikov L I 1999 Phys.-Usp. 42 107
[24] Mlynek J Abdumalikov A Eichler C Wallraff A 2014 Nat. Commun. 5 5186
[25] Jiang Y J H Jing H 2018 Chin. Phys. Lett. 35 044205
[26] Luo L F Ling Y S 1986 Chin. Phys. Lett. 3 293
[27] Kästel J Fleischhauer M 2005 Phys. Rev. 71 011804
[28] Temnov V V Woggon U 2005 Phys. Rev. Lett. 95 243602
[29] Martín-Cano D Martín-Moreno L García-Vidal F J Moreno E 2010 Nano Lett. 10 3129
[30] Choquette J Marzlin K P Sanders B 2010 Phys. Rev. 82 023827
[31] Pan J Sandhu S Huo Y Stuhrmann N Povinelli M J Harris J S Fejer M Fan S 2010 Phys. Rev. 81 041101
[32] Fleury R Alü A 2013 Phys. Rev. 87 201101
[33] Tan R Huang H 2014 Chin. Phys. Lett. 31 084205
[34] Kalachev A Kröll S 2006 Phys. Rev. 74 023814
[35] Chen Y N Chen S L Lambert N Li C M Chen G Y Nori F 2013 Phys. Rev. 88 052320
[36] Agarwal G S Kapale K T 2006 Phys. Rev. 73 022315
[37] Wang X Schirmer S G 2010 arXiv: 1005.2114v2 [quant-ph]
[38] Röhlsberger R Schlage K Sahoo B Couet S Rffer R 2010 Science 328 1248
[39] Scully M O Zubairy M S 1997 Quantum optics Cambridge Cambridge university press
[40] Agarwal G S 2012 Quantum optics Cambridge Cambridge university press
[41] Zhu C J Yang Y P Agarwal G S 2017 Phys. Rev. 95 063842